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Abstract 
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hadrons, deeply inelastic scattering, and hard processes in hadron-hadron collisions. 
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1 Introduction 



A prediction for experiment based on perturbative QCD combines a particular calculation of 
Feynman diagrams with the use of general features of the theory. The particular calculation is 
easy at leading order, not so easy at next-to-leading order and extremely difficult beyond the 
next-to-leading order. This calculation of Feynman diagrams would be a purely academic 
exercise if we did not use certain general features of the theory that allow the Feynman 
diagrams to be related to experiment: the renormalization group and the running coupling; 
the existence of infrared safe observables; the factorization property that allows us to isolate 
hadron structure in parton distribution functions. 

In these lectures, I discuss these structural features of the theory that allow a comparison 
of theory and experiment. Along the way we will discover something about certain important 
processes: e~^e~ annihilation; deeply inelastic scattering; hard processes in hadron-hadron 
collisions. By discussing the particular along with the general, I hope to arm the reader 
with information that speakers at research conferences take to be collective knowledge - 
knowledge that they assume the audience already knows. 

Now here is the disclaimer. We will not learn how to do significant calculations in QCD 
perturbation theory. Three lectures is not enough for that. 

I hope that the reader may be inspired to pursue the subjects discussed here in more 
detail. A good source is the Handbook of Perturbative QCD |T]] by the CTEQ collaboration. 
More recently, Ellis, Stirling and Webber have written an excellent book that covers the 
most of the subjects sketched in these lectures. For the reader wishing to gain a mastery 
of the theory, I can recommend the recent books on quantum field theory by Brown [Q, 
Sterman [Q, Peskin and Schroeder and Weinberg [§. Another good source, including 
both theory and phenomenology, is the lectures in the 1995 TASI proceedings, QCD and 
Beyond . I have published a substantially similar set of lectures in the proceedings of the 
1996 SLAG Summer school |]. 

2 Electron-positron annihilation and jets 

In this section, I explore the structure of the final state in QCD. I begin with the kinematics 
of e'^e~ — i> 3 partons, then examine the behavior of the cross section for e~^e~ 3 partons 
when two of the parton momenta become coUinear or one parton momentum becomes soft. 
In order to illustrate better what is going on, I introduce a theoretical tool, null-plane 
coordinates. Using this tool, I sketch a space-time picture of the singularities that we find in 
momentum space. The singularities of perturbation theory correspond to long-time physics. 
We see that the structure of the final state suggested by this picture conforms well with what 
is actually observed. 

I draw a the distinction between short-time physics, for which perturbation theory is 
useful, and long-time physics, for which the perturbative expansion is out of control. Finally, 
I discuss how certain experimental measurements can probe the short-time physics while 
avoiding sensitivity to the long-time physics. 
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2.1 Kinematics of e^e 3 partons 




Figure 1: Feynman diagram for e^e ^ qqg- 

Consider the process e"'"e~ qqg, as illustrated in Fig. |l|. Let -y/i be the total energy in 
the cm. frame and let q'^ be the virtual photon (or Z boson) momentum, so g^g^ = s. Let 
Pi be the momenta of the outgoing partons (g, q, g) and let Ei = p° be the energies of the 
outgoing partons. It is useful to define energy fractions Xi by 

Ei 2p, -q 
a/s/2 s 

< Xi. (2) 



Then 

Energy conservation gives 



Thus only two of the independent. 

Let 6ij be the angle between the momenta of partons i and j. We can relate these angles 
to the momentum fractions as follows: 

2pi ■ P2 = (Pi + P2f = (q- Pzf = s-2q-p3, (4) 

2^1^2(1 -cos ^12) =s(l-X3). (5) 

Dividing this equation by s/2 and repeating the argument for the two other pairs of partons, 
we obtain three relations for the angles 6ij: 



XiX2(l - cos6'i2) = 2(1 -X3) 

X2a;3(l - cos6'23) = 2(1 -a;i), 

X3a;i(l - cos6'3i) = 2(1-0:2). (6) 
We learn two things immediately. First, 

Xi < 1. (7) 
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Second, the three possible colhnear configurations of the partons are mapped into Xi space 
very simply: 

012^0 ^ 1, 
6*23 ^ xi 1, 

^31^0 ^ X2^l. (8) 




Figure 2: Allowed region for {xi,X2)- Then X3 is 2 — si — 2:2. The labels and small pictures in 
the right hand diagram show the physical configuration of the three partons corresponding 
to subregions in the allowed triangle. 

The relations < < 1, together with X3 = 2 — xi — X2, imply that the allowed region for 
(xi,X2) is a triangle, as shown in Fig. |^. The edges Xj = 1 of the allowed region correspond 
to two partons being collinear, as also shown in Fig. |^. The corners Xj = correspond to 
one parton momentum being soft (pf 0). 



2.2 Structure of the cross section 

One can easily calculate the cross section corresponding to Fig. [l| and the similar amplitude 
in which the gluon attaches to the antiquark line. The result is 

1_ da ^ as x\ + xl 

ao dxidx2 2n (1 — xi)(l — X2) ' 

where Cp = 4/3 and ctq = {Ana'^/s) is the total cross section for e+e" hadrons at 

order The cross section has collinear singularities: 

(1 — Xl) 0, (2&3 collinear); 

(I-X2) 0, (1&3 collinear). (10) 

There is also a singularity when the gluon is soft: X3 —>■ 0. In terms of xi and X2, this 
singularity occurs when 

(l-xi)^O, (1-X2)^0, ^^ ~^^\ r^ const. (11) 

(1 - X2j 
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Let us write the cross section in a way that displays the coUinear singularity at ^31 
and the soft singularity at E3 ^ 0: 



1 



da 



(To dE^ d cos ^^31 



2tt ^ ^3(1 - cos^si) 



(12) 



Here f{E^, 631) a rather complicated function. The only thing that we need to know about 
it is that it is finite for E^ —>■ and for 6'3i 0. 

Now look at the collinear singularity, ^31 ^0. If we integrate over the singular region 
holding E^ fixed we find that the integral is divergent: 



d cos 6*31 



da 



dE^ d cos ^31 



log(oo). 



(13) 



Similarly, if we integrate over the region of the soft singularity, holding ^31 fixed, we find 
that the integral is divergent: 



dEs = log(oo). 

dE^ d cos &3i 



(14) 



Evidently, perturbation theory is telling us that we should not take the perturbative cross 
section too literally. The total cross section for e~^e~ hadrons is certainly finite, so this 
partial cross section cannot be infinite. What we are seeing is a breakdown of perturbation 
theory in the soft and collinear regions, and we should understand why. 




Figure 3: Cross section for e^e 
or collinear with the quark. 



qqg., illustrating the singularity when the gluon is soft 



Where do the singularities come from? Look at Fig. ^ (in a physical gauge). The 
scattering matrix element M. contains a factor l/(pi +^3)^ where 

{Pi+Pz? = 2pi ■ P3 = 2^1^3(1 - cos^3i)- (15) 

Evidently, l/(pi +^3)^ is singular when 6^1 — > and when E3 0. The collinear singularity 
is somewhat softened because the numerator of the Feynman diagram contains a factor 
proportional to 631 in the collinear limit. (This is not exactly obvious, but is easily seen by 
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calculating. If you like symmetry arguments, you can derive this factor from quark helicity 
conservation and overall angular momentum conservation.) We thus find that 



\M\^ oc 
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1 2 



31 



3(731 



(16) 



for E3 ^ and ^31 — 0. Note the universal nature of these factors. 

Integration over the double singular region of the momentum space for the gluon has the 
form 



E^dE^d cos 9si( 
E-i 



E-^dE-^de. 



31^ 



(17) 



Combining the integration with the matrix element squared gives 

1 2 

JTT^ J/l2 J J, f^31 

^3 



da 



E^dEsdOl^c 



dE^ '^^31 



E, 9^ 



31 



Thus we have a double logarithmic divergence in perturbation theory for the soft and coUinear 
region. With just a little enhancement of the argument, we see that there is a collinear diver- 
gence from integration over ^31 at finite E^ and a separate soft divergence from integration 
over E3 at finite ^31. Essentially the same argument applies to more complicated graphs. 
There are divergences when two final state partons become collinear and when a final state 
gluon becomes soft. Generalizing further 0, there are also divergences when several final 
state partons become collinear to one another or when several (with no net flavor quantum 
numbers) become soft. 

We have seen that if we integrate over the singular region in momentum space with no 
cutoff, we get infinity. The integrals are logarithmically divergent, so if we integrate with 
an infrared cutoff Mj^, we will get big logarithms of Mf^/s. Thus the collinear and soft 
singularities represent perturbation theory out of control. Carrying on to higher orders of 
perturbation theory, one gets 



1 + X (big) + X (big)^ + ■ ■ 



(19) 



If this expansion is in powers of as{Mz), we have <C 1. Nevertheless, the big logarithms 
seem to spoil any chance of the low order terms of perturbation theory being a good ap- 
proximation to any cross section of interest. Is the situation hopeless? We shall have to 
investigate further to see. 



2.3 Interlude: Null plane coordinates 

In order to understand better the issue of singularities, it is helpful to introduce a concept that 
is generally quite useful in high energy quantum field theory, null plane coordinates. The idea 
is to describe the momentum of a particle using momentum components = {p'^ , p~ , , p^) 
where 

p± = (/±p=^)/v^. (20) 
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pO 




Figure 4: Null plane axes in momentum space. 



For a particle with large momentum in the +z direction and limited transverse momentum, 
p"*" is large and p~ is small. Often one chooses the plus axis so that a particle or group of 
particles of interest have large and small p~ and pt- 

Using null plane components, the covariant square of p^ is 

p2 = 2p+p" - p^. (21) 
Thus, for a particle on its mass shell, p~ is 

(-) 

Note also that, for a particle on its mass shell, 

j9+ > , > . (23) 

Integration over the mass shell is 

<-)-V5;|fe.-(-)-v<'vf|;.... (.4) 

We also use the plus/minus components to describe a space-time point x^: = (x°± 
x^)/^/2. In describing a system of particles moving with large momentum in the plus direc- 
tion, we are invited to think of x~^ as "time." Classically, the particles in our system follow 
paths nearly parallel to the x~^ axis, evolving slowly as it moves from one x~^ = const, plane 
to another. 

We relate momentum space to position space for a quantum system by Fourier trans- 
forming. In doing so, we have a factor exp{ip ■ x), which has the form 

p ■ X = p'^x~ + p~x^ — Pt ■ x^. (25) 

Thus x~ is conjugate to p~^ and x"*" is conjugate to p~ . That is a little confusing, but it is 
simple enough. 
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Figure 5: Correspondence between singularities in momentum space and the development 
of the system in space-time. 



2.4 Space-time picture of the singularities 

We now return to the singularity structure of e~^e^ qqg. Define pi + = k^- Choose 
null plane coordinates with k~^ large and k^- = 0. Then /c^ = 2k^k^ becomes small when 

2pt 2pt ^ ^ 

becomes small. This happens when ps^^ becomes small with fixed pf and p^ , so that the 
gluon momentum is nearly collinear with the quark momentum. It also happens when p3 ^ 
and pf both become small with p^ cx |p3,t|, so that the gluon momentum is soft. ( It also 
happens when the quark becomes soft, but there is a numerator factor that cancels the soft 
quark singularity.) Thus the singularities for a soft or collinear gluon correspond to small 
k-. 

Now consider the Fourier transform to coordinate space. The quark propagator in Fig. ^ 

is r _ 

Splk) = / dx^dx d:K exp{i[k~^ X +k x~^ — k ■ x\) Sf{x). (27) 

When k~^ is large and k~ is small, the contributing values of x have small x~ and large x~^. 
Thus the propagation of the virtual quark can be pictured in space-time as in Fig. |^. The 
quark propagates a long distance in the x+ direction before decaying into a quark-gluon pair. 
That is, the singularities that can lead to divergent perturbative cross sections arise from 
interactions that happen a long time after the creation of the initial quark- ant iquark pair. 



2.5 Nature of the long-time physics 

Imagine dividing the contributions to a scattering cross section into long-time contributions 
and short-time contributions. In the long-time contributions, perturbation theory is out 
of control, as indicated in Eq. ([T9|). Nevertheless the generic structure of the long-time 
contribution is of great interest. This structure is illustrated in Fig. ^. Perturbative diagrams 
have big contributions from space-time histories in which partons move in collinear groups 
and additional partons are soft and communicate over large distances, while carrying small 
momentum. 



8 



Figure 6: Typical paths of partons in space contributing to e~^e~ — > hadrons, as suggested 
by the singularities of perturbative diagrams. Short wavelength fields are represented by 
classical paths of particles. Long wavelength fields are represented by wavy lines. 



The picture of Fig. ^ is suggested by the singularity structure of diagrams at any fixed 
order of perturbation theory. Of course, there could be nonperturbative effects that would 
invalidate the picture. Since nonperturbative effects can be invisible in perturbation theory, 
one cannot claim that the structure of the final state indicated in Fig. |^ is known to be a 
consequence of QCD. One can point, however, to some cases in which one can go beyond 
fixed order perturbation theory and sum the most important effects of diagrams of all orders 
(for example, Ref. [|10|)- In such cases, the general picture suggested by Fig. | remains intact. 



We thus find that perturbative QCD suggests a certain structure of the final state pro- 
duced in e~^e~ — > hadrons: the final state should consist of jets of nearly coUinear particles 
plus soft particles moving in random directions. In fact, this qualitative prediction is a 
qualitative success. 

Given some degree of qualitative success, we may be bolder and ask whether perturbative 
QCD permits quantitative predictions. If we want quantitative predictions, we will somehow 
have to find things to measure that are not sensitive to interactions that happen long after 
the basic hard interaction. This is the subject of the next section. 



2.6 The long-time problem 

We have seen that perturbation theory is not effective for long-time physics. But the detector 
is a long distance away from the interaction, so it would seem that long-time physics has to 
be present. 

Fortunately, there are some measurements that are not sensitive to long-time physics. 
An example is the total cross section to produce hadrons in e~^e~ annihilation. Here effects 
from times At S> I/a/s cancel because of unitarity. To see why, note that the quark state 
is created from the vacuum by a current operator J at some time t; it then develops from 
time t to time oo according to the interaction picture evolution operator U{oo,t), when it 
becomes the final state |A^). The cross section is proportional to the sum over of this 
amplitude times a similar complex conjugate amplitude with t replaced by a different time 
t'. We Fourier transform this with exp{—i^/s {t — t')), so that we can take At = t — t' to be 
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of order Now replacing J2 by the unit operator and using the unitarity of the 

evolution operators U, we obtain 



^(0|J(Ot/(t',oo)|iV)(iV|f/(oo,t)J(t)|0) 



(28) 



N 



= {0\J{t')U{t',oo)U{oo,t)J{t)\0) = {0\J{t')U{t',t)J{t)\0). 



Because of unitarity, the long-time evolution has canceled out of the cross section, and we 
have only evolution from t to t'. 

There are three ways to view this result. First, we have the formal argument given 
above. Second, we have the intuitive understanding that after the initial quarks and gluons 
are created in a time At of order something will happen with probability 1. Exactly 

what happens is long-time physics, but we don't care about it since we sum over all the 
possibilities |A^). Third, we can calculate at some finite order of perturbation theory. Then 
we see infrared infinities at various stages of the calculations, but we find that the infinities 
cancel between real gluon emission graphs and virtual gluon graphs. An example is shown 
in Fig. 0. 




Figure 7: Cancellation between real and virtual gluon graphs. If we integrate the real gluon 
graph on the left times the complex conjugate of the similar graph with the gluon attached 
to the antiquark, we will get an infrared infinity. However the virtual gluon graph on the 
right times the complex conjugate of the Born graph is also divergent, as is the Born graph 
times the complex conjugate of the virtual gluon graph. Adding everything together, the 
infrared infinities cancel. 

We see that the total cross section is free of sensitivity to long-time physics. If the total 
cross section were all you could look at, QCD physics would be a little boring. Fortunately, 
there are other quantities that are not sensitive to infrared effects. They are called infrared 
safe quantities. 

To formulate the concept of infrared safety, consider a measured quantity that is con- 
structed from the cross sections. 



to make n hadrons in e'^e annihilation. Here Ej is the energy of the jth hadron and 
flj = {9j, describes its direction. We treat the hadrons as effectively massless and do not 



(icr[n] 



(29) 



dVt2dE'idVt'i ■ ■ ■ dEndVt, 



•n 
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distinguish the hadron flavors. Following the notation of Ref. let us specify functions 
Sn that describe the measurement we want, so that the measured quantity is 

+ ^ J dVL2dE^dQ.'idE^dVLi 



dn^dEsdn^dE.dn, ^'^p'^p'^^p'^ 



+ •••• (30) 

The functions S are symmetric functions of their arguments. In order for our measurement 
to be infrared safe, we need 

Sn+M, . . . , (1 - AK, AK) = Snipl . . . (31) 

for < A < 1. 





Figure 8: Infrared safety. In an infrared safe measurement, the three jet event shown on the 
left should be (approximately) equivalent to an ideal three jet event shown on the right. 



What does this mean? The physical meaning is that the functions Sn and Sn-i are 
related in such a way that the cross section is not sensitive to whether or not a mother 
particle divides into two collinear daughter particles that share its momentum. The cross 
section is also not sensitive to whether or not a mother particle decays to a daughter particle 
carrying all of its momentum and a soft daughter particle carrying no momentum. The 
cross section is also not sensitive to whether or not two collinear particles combine, or a 
soft particle is absorbed by a fast particle. All of these decay and recombination processes 
can happen with large probability in the final state long after the hard interaction. But, 
by construction, they don't matter as long as the sum of the probabilities for something to 
happen or not to happen is one. 

Another version of the physical meaning is that for an IR-safe quantity a physical event 
with hadron jets should give approximately the same measurement as a parton event with 
each jet replaced by a parton, as illustrated in Fig. ||. To see this, we simply have to delete 
soft particles and combine collinear particles until three jets have become three particles. 
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In a calculation of the measured quantity X, we simply calculate with partons instead of 
hadrons in the final state. The calculational meaning of the infrared safety condition is that 
the infrared infinities cancel. The argument is that the infinities arise from soft and collinear 
configurations of the partons, that these configurations involve long times, and that the time 
evolution operator is unitary. 

I have started with an abstract formulation of infrared safety. It would be good to have 
a few examples. The easiest is the total cross section, for which 

5„K,...,p^:) = l. (32) 

A less trivial example is the thrust distribution. One defines the thrust 7^ of an particle 
event as 

r„(rf, . . . , pt:) = max . (33) 

Here m is a unit vector, which we vary to maximize the sum of the absolute values of the 
projections of pi on u. Then the thrust distribution (l/atot) da/dT is defined by taking 

SM, . . . ,P^:) = {l/crtot) 6{T - TM. ■ ■ ■,P'n)) ■ (34) 

It is a simple exercise to show that the thrust of an event is not affected by collinear parton 
splitting or by zero momentum partons. Therefore the thrust distribution is infrared safe. 
Another example is the energy-energy correlation function dT,/dcos{6) [0: 

'5n(p^ . . . ,pt:) = E ^ (cos(%) - <^os{e)) . (35) 

ij ^ 

This measures the correlation between the energies measured by detectors separated by an 
angle 9 as depicted in Fig. ^ Is this infrared safe? Note that the contribution from a particle 
with Ei ^ drops out. In addition, replacing one particle by two collinear particles doesn't 
change the thrust: 

il~X)EnEj + XEn Ej =EnEj. (36) 

This works for the autocorrelation term too: 

(1 - X)' El + 2A(1 -X)El + X' El = El (37) 

A final example is the cross section to make n jets, cr„. Intuitively, a jet is supposed to be 
a spray of particles all going in approximately the same direction. To make this precise, we 
need a definite algorithm. There are several algorithms to choose from. Here is the simplest 
(but not the best) one. 

Start with a list of momenta Pi,P2, ■ ■ ■ ,Pn- the start, these represent the momenta 
of particles. (In a perturbative calculation, they are the momenta of partons.) Choose a 
parameter j/cut- Now proceed through the following steps: 

1. Find the pair such that {pi + Pj)'^ is the smallest. 
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Figure 9: The energy-energy correlation function 





Figure 10: Jet definition 

2. If {pi + PjY > VcntS, exit. Else continue. 

3. Replace the two momenta pi and pj in the list by their sum Pk=Pi+ Pj- 

4. Go to 1. 

This produces a list of momenta pi of jets, an is the cross section to have n jets. The infrared 
safety of (T„ is easy to prove given our experience with the previous examples. 

There are several variations on this theme. For instance, change the resolution condition 
or the combination prescription. A comparison of methods can be found in Ref. |T^. 1 



discuss jet cross sections for hadron collisions in Sec. p^ . 

Before leaving this subject, 1 should mention another way to eliminate sensitivity to 
long-time physics. Consider the cross section 



(38) 



This cross section can be written as a convolution of two factors. The first factor is a 
calculated "hard scattering cross section" for e~^e~ — quark + X or e~^e~ —>■ gluon + X . The 
second factor is a "parton decay function" for quark ^ vr + X or gluon — > tt + X. These 
functions contain the long-time sensitivity and are to be measured, since they cannot be 
calculated perturbatively. However, once they are measured in one process, they can be used 
for another process. This final state factorization is similar to the initial state factorization 
involving parton distribution functions, which we will discuss later. (See Refs. |llIl,[H,|l3 ^'^'^ 
more information.) 
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3 The smallest time scales 



In this section, I explore the physics of time scales smaller than l/i/i- One way of looking at 
this physics is to say that it is plagued by infinities and we can manage to hide the infinities. 
A better view is that the short-time physics contains wonderful truths that we would like 
to discover - truths about grand unified theories, quantum gravity and the like. However, 
quantum field theory is arranged so as to effectively hide the truth from our experimental 
apparatus, which can probe with a time resolution of only an inverse half TeV. 

I first outline what renormalization does to hide the ugly infinities or the beautiful truth. 
Then I describe how renormalization leads to the running coupling. Because of renormaliza- 
tion, calculated quantities depend on a renormalization scale. I look at how this dependence 
works and how the scale can be chosen. Finally, I discuss how one can use experiment to 
look for the hidden physics beyond the Standard Model, taking high Et jet production in 
hadron collisions as an example. 



3.1 What renormalization does 

In any Feynman graph, one can insert perturbative corrections to the vertices and the prop- 
agation of particles, as illustrated in Fig. 0. The loop integrals in these graphs will get big 
contributions from momenta much larger than ^/s. That is, there are big contributions from 
interactions that happen on time scales much smaller than 1 / ^/s. 1 have tried to illustrate 
this in the figure. The virtual vector boson propagates for a time l/y/s, while the virtual 
fluctuations that correct the electroweak vertex and the quark propagator occur over a time 
At that can be much smaller than 1 / y/s. 

Let us pick an ultraviolet cutoff M that is much larger than ^/s, so that we calculate the 
effect of fluctuations with 1/M < At exactly, up to some order of perturbation theory. What, 
then, is the effect of virtual fluctuations on smaller time scales. At with At < 1/M but, say. 
At still larger than tpiank, where gravity takes over? Let us suppose that we are willing to 
neglect contributions to the cross section that are of order y/s/M or smaller compared to the 
cross section itself. Then there is a remarkable theorem |15| : the effects of the fluctuations are 
not particularly small, but they can be absorbed into changes in the couplings of the theory. 
(There are also changes in the masses of the theory and adjustments to the normalizations 
of the field operators, but we can concentrate on the effect on the couplings.) 

The program of absorbing very short-time physics into a few parameters goes under the 
name of renormalization. There are several schemes available for renormalizing. Each of 
them involves the introduction of some scale parameter that is not intrinsic to the theory 
but tells how we did the renormalization. Let us agree to use MS renormalization (see 
Ref. for details). Then we introduce an MS renormalization scale /i. A good (but 



approximate) way of thinking of /i is that the physics of time scales At <^ is removed 
from the perturbative calculation. The effect of the small time physics is accounted for by 
adjusting the value of the strong coupling, so that its value depends on the scale that we 
used: = a<((/i). (The value of the electromagnetic coupling also depends on fi.) 
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Figure 11: Renormalization. The effect of the very small time interactions pictured are 
absorbed into the running coupling. 



3.2 The running coupling 



Tnnnnni(J)innnnnr 




Figure 12: Short-time fluctuations in the propagation of the gluon field absorbed into the 
running strong coupling. 

We account for time scales much smaller than 1/yU by using the running coupling as(/i). 
That is, a fluctuation such as that illustrated in Fig. can be dropped from a calculation and 
absorbed into the running coupling that describes the probability for the quark in the figure 
to emit the gluon. The /i dependence of Q;s(yu) is given by a certain differential equation, 
called the renormalization group equation (see Ref. p!5[]): 



pM,)) = -0j^Y-0j^U---. (39) 



(iln(/i2) TT \ ^ 



TT 



One calculates the beta function /?(«<() perturbatively in QCD. The first coefficient, with the 
conventions used here, is 

/?o = (33-2iV^)/12, (40) 

where Nf is the number of quark flavors. 

Of course, at time scales smaller than a very small cutoff 1/M (at the "GUT scale," say) 
there is completely different physics operating. Therefore, if we use just QCD to adjust the 
strong coupling, we can say that we are accounting for the physics between times 1/M and 
l//i. The value of at /iq ~ M is then the boundary condition for the differential equation. 



See Fig. |3 



The renormalization group equation sums the effects of short-time fluctuations of the 
fields. To see what one means by "sums" here, consider the result of solving the renormal- 
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renomalization 

(jrouD fixed order 
1 



log(l/M) log(l/n) 



log{At) 



Figure 13: Time scales accounted for by explicit fixed order perturbative calculation and by 
use of the renormalization group. 

ization group equation with all of the jSi beyond /5o set to zero: 

\2 1^2/2 //i#2n ^3/ 



+ (/?o/vr)'ln^(/x7^ ) ai{M) + 

as{M) 

l + (/3o/7r) a,(M) ln(/i7M2)- 



(41) 



A series in powers of as{M) - that is the strong coupling at the GUT scale - is summed into 
a simple function of fi. Here as{M) appears as a parameter in the solution. 

Note a crucial and wonderful fact. The value of ^^(/i) decreases as increases. This 
is called "asymptotic freedom." Asymptotic freedom implies that QCD acts like a weakly 
interacting theory on short time scales. It is true that quarks and gluons are strongly bound 
inside nucleons, but this strong binding is the result of weak forces acting collectively over a 
long time. 

In Eq. (^TD, we are invited to think of the graph of as(/i) versus fi. The differential equa- 
tion that determines this graph is characteristic of QCD. There could, however, be different 
versions of QCD with the same differential equation but different curves, corresponding to 
different boundary values ^^(M). Thus the parameter as{M) tells us which version of QCD 
we have. To determine this parameter, we consult experiment. Actually, Eq. (^) is not the 
most convenient way to write the solution for the running coupling. A better expression is 

Here we have replaced as{M) by a different (but completely equivalent) parameter A. A 
third form of the running coupling is 

Here the value of as{fi) at /i = Mz labels the version of QCD that obtains in our world. 

In any of the three forms of the running coupling, one should revise the equations to 
account for the second term in the beta function in order to be numerically precise. 

3.3 The choice of scale 

In this section, we consider the choice of the renormalization scale /i in a calculated cross 
section. Consider, as an example, the cross section for e~^e~ — >■ hadrons via virtual photon 
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decay. Let us write this cross section in the form 

^tot = ^(|EQ?)[l + A]. (44) 

Here s is the square of the cm. energy, a is e^/(47r), and Qj is the electric charge in units of 
e carried by the quark of flavor /, with f = u,d, s, c, b. The nontrivial part of the calculated 
cross section is the quantity A, which contains the effects of the strong interactions. Using 
MS renormalization with scale fi, one finds (after a lot of work) that A is given by Ref. ||16||: 



' [1.4092 + 1.9167 In (/.V^)! 



+ 



TT 



TT 



-12.805 + 7.8186 In (^V^) + 3.674 ln2(^Vs^^ 

+ ■■■. ' (45) 

Here, of course, one should use for a;s(/i) the solution of the renormalization group equation 
(^) with at least two terms included. 

As discussed in the preceding subsection, when we renormalize with scale /i, we are 
defining what we mean by the strong coupling. Thus as in Eq. (^) depends on /i. The 
perturbative coefficients in Eq. ( |45| ) also depend on /i. On the other hand, the physical cross 
section does not depend on fx: 

^ = 0- (46) 

a m jjL'^ 

That is because /x is just an artifact of how we organize perturbation theory, not a parameter 
of the underlying theory. 

Let us consider Eq. (^6[) in more detail. Write A in the form 



A ~ E c„(/i) «.(/i)". (47) 

n=l 

If we differentiate not the complete infinite sum but just the first N terms, we get minus the 
derivative of the sum from + 1 to infinity. This remainder is of order a^^^ as as — >■ 0. 
Thus 

d ^ 

——^ Y: cn(/i) as{^^r ~ 0{as{^^r^'). (48) 

"^^A'' n=l 

That is, the harder we work calculating more terms, the less the calculated cross section 
depends on /i. 

Since we have not worked infinitely hard, the calculated cross section depends on /i. What 
choice shall we make for /i? Clearly, In (/i^/ s) should not be big. Otherwise the coefficients 
Cn(/x) are large and the "convergence" of perturbation theory will be spoiled. There are 
some who will argue that one scheme or the other for choosing /i is the "best." You are 
welcome to follow whichever advisor you want. I will show you below that for a well behaved 
quantity like A the precise choice makes little difference, as long as you obey the common 
sense prescription that ln(/i^/s) not be big. 
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3.4 An example 
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Figure 14: Dependence of A(/i) on the MS renormalization scale /z. The falling curve is Ai. 
The flatter curve is A2. The horizontal hues indicates the amount of variation of A2 when 
fi varies by a factor 2. 



Let us consider a quantitative example of how A(yu) depends on fi. This will also give us 
a chance to think about the theoretical error caused by replacing A by the sum A„ of the 
flrst n terms in its perturbative expansion. Of course, we do not know what this error is. 
All we can do is provide an estimate. (Our discussion will be rather primitive. For a more 



detailed error estimate for the case of the hadronic width of the Z boson, see Ref. |T7| .) 

Let us think of the error estimate in the spirit of a "1 a" theoretical error: we would be 
surprised if | A„ — A | were much less than the error estimate and we would also be surprised 
if this quantity were much more than the error estimate. Here, one should exercise a little 
caution. We have no reason to expect that theory errors are gaussian distributed. Thus 
a 4 (J difference between A„ and A is not out of the question, while a 4ct fluctuation in a 
measured quantity with purely statistical, gaussian errors is out of the question. 

Take as(M^) = 0.117, ^/s = 34 GeV, 5 flavors. In Fig. |1^, 1 plot A(yu) versus p deflned 

by 

/i = (49) 

The steeply falling curve is the order a] approximation to A(/i), Ai(yu) = as{f^)/n. Notice 
that if we change /i by a factor 2, Ai(/i) changes by about 0.006. If we had no other 
information than this, we might pick Ai(y^) ^ 0.044 as the "best" value and assign a 
±0.006 error to this value. (There is no special magic to the use of a factor of 2 here. The 
reader can pick any factor that seems reasonable.) 

Another error estimate can be based on the simple expectation that the coefficients of 
a" are of order 1 for the flrst few terms. (Eventually, they will grow like n\. Ref. Jl^ takes 
this into account, but we ignore it here.) Then the flrst omitted term should be of order 
±1 X ctg ^ ±0.020 using as(34 GeV) ~ 0.14. Since this is bigger than the previous ±0.006 
error estimate, we keep this larger estimate: A ^ 0.044 ± 0.020. 
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Returning now to Fig. the second curve is the order approximation, A2(/u). Note 
that A2(/i) is less dependent on fi than Ai(yu). 

What value would we now take as our best estimate of A? One idea is to choose the 
value of fi at which A2(/i) is least sensitive to /i. This idea is called the principle of minimal 
sensitivity p!8| : 



A 



PMS 



A(/i 



PMS) 



dln/i 



0. 



(50) 



This prescription gives A ^ 0.0470. Note that this is about 0.003 away from our previous 
estimate, A ^ 0.0440. Thus our previous error estimate of 0.020 was too big, and we should 
be surprised that the result changed so little. We can make a new error estimate by noting 
that A2(/i) varies by about 0.0012 when /i changes by a factor 2 from Upms- Thus we might 
estimate that A ^ 0.0470 with an error of ±0.0012. This estimate is represented by the two 
horizontal lines in Fig. IT^. 

An alternative error estimate can be based on the next term being of order ±1 x 
a^(34 GeV) ~ 0.003. Since this is bigger than the previous ±0.0012 error estimate, we 
keep this larger estimate: A ^ 0.0470 ± 0.003. 

I should emphasize that there are other ways to pick the "best" value for A. For instance, 
one can use the BLM method [|1^], which is based on choosing the fi that sets to zero the 
coefficient of the number of quark flavors in A2(/i). Since the graph of A2(/u) is quite flat, it 
makes very little difference which method one uses. 




Figure 15: Dependence of A(/i) on the MS renormalization scale fi, flrst with a normal scale 
and then with an expanded scale. The falling curve is Ai. The flatter curve is A2. The still 
flatter curve is A3. The horizontal lines represent the variation of A2 when /i varies by a 
factor 2. 

Now let us look at A(/i) evaluated at order a^, A^lfi). Here we make use of the full 
formula in Eq. (^). In Fig. |T^, I plot A^^^) along with A2(/i) and Ai(/i). The variation of 
A3(/i) with /i is smaller than that of A2(/i). The improvement is not overwhelming, but is 
apparent particularly at small /x. 
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It is a little difficult to see what is happening in the ffist graph of Fig. |T^, so I show the 
same thing with an expanded scale. (Here the error band based on the /i dependence of A2 
is also indicated. Recall that we decided that this error band was an underestimate.) The 
curve for A^lfi) has zero derivative at two places. The corresponding values are A ^ 0.0436 
and A f» 0.0456. If I take the best value of A to be the average of these two values and the 
error to be half the difference, I get A ^ 0.0446 ± 0.0010. 

The alternative error estimate is ±1 x a^(34 GeV) ^ 0.0004. We keep the lar ger error 
estimate of ±0.0010. 

Was the previous error estimate valid? We guessed A ^ 0.0470 ± 0.003. Our new 
best estimate is 0.0446. The difference is 0.0024, which is in line with our previous error 
estimate. Had we used the error estimate ±0.0012 based on the fi dependence, we would 
have underestimated the difference, although we would not have been too far off. 



3.5 Beyond the Standard Model 

We have seen how the renormalization group enables us to account for QCD physics at time 



scales much smaller than ^/s, as indicated in Fig. 0. However, at some scale At ~ 1/ M, we 
run into the unknown! 

How can we see the unknown in current experiments? First, the unknown physics affects 
as, aem, sin^(6'iy)- Second, the unknown physics affects masses of u,d, . . . ,e, fi, . . .. That 
is, the unknown physics (presumably) determines the parameters of the Standard Model. 
These parameters have been well measured. Thus, a Nobel prize awaits the physicist who 
figures out how to use a model for the unknown physics to predict these parameters. 





Figure 16: New physics at a TeV scale. In the ffist diagram, quarks scatter by gluon 
exchange. In the second diagram, the quarks exchange a new object with a TeV mass, or 
perhaps exchange some of the constituents out of which quarks are made. 



There is another way that as yet unknown physics can affect current experiments. Sup- 
pose that quarks can scatter by the exchange of some new particle with a heavy mass M, 
as illustrated in Fig. 0, and suppose that this mass is not too enormous, only a few TeV. 
Perhaps the new particle isn't a particle at all, but is a pair of constituents that live inside 
of quarks. As mentioned above, this physics affects the parameters of the Standard Model. 
However, unless we can predict the parameters of the Standard Model, this effect does not 
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help us. There is, however, another possible clue. The physics at the TeV scale can introduce 
new terms into the lagrangian that we can investigate in current experiments. 



In the second diagram in Fig. |T6|, the two vertices are never at a separation in time 
greater than 1/M, so that our low energy probes cannot resolve the details of the structure. 
As long as we stick to low energy probes, ^/s -C M, the effect of the new physics can be 
summarized by adding new terms to the lagrangian of QCD. A typical term might be 

A/: = ^ ^ri^ ^TmV'- (51) 

There is a factor that represents how well the new physics couples to quarks. The most 
important factor is the factor 1/M^. This factor must be there: the product of field operators 
has dimension 6 and the lagrangian has dimension 4, so there must be a factor with dimension 
—2. Taking this argument one step further, the product of field operators in AC must have 
a dimension greater than 4 because any product of field operators having dimension equal 
to or less than 4 that respects the symmetries of the Standard Model is already included in 
the lagrangian of the Standard Model. 

3.6 Looking for new terms in the effective lagrangian 

How can one detect the presence in the lagrangian of a term like that in Eq. (|5T|) ? These 
terms are small. Therefore we need either a high precision experiment, or an experiment 
that looks for some effect that is forbidden in the Standard Model, or an experiment that 
has moderate precision and operates at energies that high as possible. 

Let us consider an example of the last of these possibilities, p + p — >• jet + X as a function 
of the transverse energy (~ Pt) of the jet. The new term in the lagrangian should add a 
little bit to the observed cross section that is not included in the standard QCD theory. 
When the transverse energy Et of the jet is small compared to M, we expect 

Data — Theory _2 -^r ^^2) 
Theory 

Here the factor g'^/M'^ follows because AC contains this factor. The factor follows because 
the left hand side is dimensionless and Et is the only factor with dimension of mass that is 
available. 

In Fig. |l^, I show a plot comparing experimental jet cross sections from CDF and 



DO 1^ compared to next-to-leading order QCD theory. The theory works fine for Et < 
200 GeV, but for 200 GeV < Et, there appears to be a systematic deviation of just the form 
anticipated in Eq. 

This example illustrates the idea of how small distance physics beyond the Standard 
Model can leave a trace in the form of small additional terms in the effective lagrangian 
that controls physics at currently available energies. However, in this case, there is some 
indication that the observed effect might be explained by some combination of the experi- 
mental systematic error and the uncertainties inherent in the theoretical prediction [^]. In 
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Figure 17: Jet cross sections from CDF and DO compared to QCD theory. (Data — The- 
ory)/Theory is plotted versus the transverse energy of the jet. The theory here is next- 
to-leading order QCD using the CTEQ3M parton distribution. Source: Ref. pO| 



particular, the prediction is sensitive to the distributions of quarks and gluons contained in 
the colliding protons, and the gluon distribution in the kinematic range of interest here is 
rather poorly known. In the next section, we turn to the definition, use, and measurement 
of the distributions of quarks and gluons in hadrons. 



4 Deeply inelastic scattering 

Until now, I have concentrated on hard scattering processes with leptons in the initial state. 
For such processes, we have seen that the hard part of the process can be described using 
perturbation theory because ^^(/i) gets small as /i gets large. Furthermore, we have seen how 
to isolate the hard part of the interaction by choosing an infrared safe observable. But what 
about hard processes in which there are hadrons in the initial state? Since the fundamental 
hard interactions involve quarks and gluons, the theoretical description necessarily involves 
a description of how the quarks and gluons are distributed in a hadron. Unfortunately, 
the distribution of quarks and gluons in a hadron is controlled by long-time physics. We 
cannot calculate the relevant distribution functions perturbatively (although a calculation 
in lattice QCD might give them, in principle). Thus we must find how to separate the 
short-time physics from the parton distribution functions and we must learn how the parton 
distribution functions can be determined from the experimental measurements. 

In this section, I discuss parton distribution functions and their role in deeply inelastic 
lepton scattering (DIS). This includes e+p — > e + X and i>+p — > e + X where the momentum 
transfer from the lepton is large. I first outline the kinematics of deeply inelastic scattering 
and define the structure functions Fi, F2 and F3 used to describe the process. By examining 
the space-time structure of DIS, we will see how the cross section can be written as a 
convolution of two factors, one of which is the parton distribution functions and the other 
of which is a cross section for the lepton to scatter from a quark or gluon. This factorization 
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involves a scale jj^p that, roughly speaking, divides the soft from the hard regime; I discuss 
the dependence of the calculated cross section on fip. With this groundwork laid, I give the 
MS definition of parton distribution functions in terms of field operators and discuss the 
evolution equation for the parton distributions. I close the section with some comments on 
how the parton distributions are, in practice, determined from experiment. 



4.1 Kinematics of deeply inelastic lepton scattering 




In deeply inelastic scattering, a lepton with momentum fc^ scatters on a hadron with 
momentum p^^. In the final state, one observes the scattered lepton with momentum k'^ as 
illustrated in Fig. 0. The momentum transfer 

q^' = k^'- k'^" (53) 

is carried on a photon, or aW or Z boson. 

The interaction between the vector boson and the hadron depends on the variables g'^ 
and p'^. From these two vectors we can build two scalars (not counting = p^). The first 
variable is 

= -q\ (54) 

where the minus sign is included so that is positive. The second scalar is the dimensionless 
Bjorken variable, 

Xhi = • (55) 

2p ■ q 

(In the case of scattering from a nucleus containing A nucleons, one replaces p^ by p^^/A and 
defines Xbj = AQV(2p-g).) 

One calls the scattering deeply inelastic if is large compared to 1 GeV^. Traditionally, 
one speaks of the scaling limit, Q"^ ^ oo with Xbj fixed. Actually, the asymptotic theory to 
be described below works pretty well if is bigger than, say, 4 GeV^ and Xbj is anywhere 
in the experimentally accessible range, roughly 10~^ < Xbj < 0.5. 
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The invariant mass squared of the hadronic final state is W"^ = {p + q)"^- In the scahng 
regime of large one has 

This justifies saying that the scattering is not only inelastic but deeply inelastic. 

We have spoken of the scalar variables that one can form from p'^ and q^^. Using the 
lepton momentum k^^, one can also form the dimensionless variable 

4.2 Structure functions for DIS 

One can make quite a lot of progress in understanding the theory of deeply inelastic scattering 
without knowing anything about QCD except its symmetries. One expresses the cross section 
in terms of three structure functions, which are functions of Xbj and only. 

Suppose that the initial lepton is a neutrino, t'^, and the final lepton is a muon. Then in 
Fig. |18| the exchanged vector boson, call it V , is a,W boson, with mass My = Mw- Alter- 
natively, suppose that both the initial and final leptons are electrons and let the exchanged 
vector boson be a photon, with mass My = 0. This was the situation in the original DIS 
experiments at SLAG in the late 1960 's. In experiments with sufficiently large Q^, Z boson 
exchange should be considered along with photon exchange, and the formalism described 
below must be augmented. 

Given only the electroweak theory to tell us how the vector boson couples to the lepton, 
one can write the cross section in the form 

where Cy is 1 in the case that is a photon and 1/(64 sin^^v^) in the case that is a 
boson. The tensor L'^'^ describes the lepton coupling to the vector boson and has the form 

= Iti {k ■ 7 -f^k' ■ 7 (59) 

in the case that is a photon. For a W boson, one has 

L^"" = Tr (A; • 7 Vk' ■ T") , (60) 

where is 7''(l-75) for a boson {u W^£) or 7^(1+75) for a W~ boson (z7 W^i). 
See Ref. [§. 

The tensor W^" describes the coupling of the vector boson to the hadronic system. It 
depends on and q^. We know that it is a Lorentz tensor and that W^^ = W^"*. We 
also know that the current to which the vector boson couples is conserved (or in the case of 
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the axial current, conserved in the absence of quark masses, which we here neglect) so that 
q^W^^ = 0. Using these properties, one finds three possible tensor structures for W^'^. Each 
of the three tensors multiplies a structure function, Fi, F2 or F3, which, since it is a Lorentz 
scalar, can depend only on the invariants Xbj and Q^. Thus 



W, 



p ■ g\ 1 
V J p-q 



F2(xbj,g') 



-i^tiuXaP^q'' F3(Xbj,(5^). 

p-q 



(61) 



If we combine Eqs. ( pq , |59| , |60| , |6Tl) , we can write the cross section for deeply inelastic 
scattering in terms of the three structure functions. Neglecting the hadron mass compared 
to Q^, the result is 



da 



(ia;bj dy 

Here the normalization factor N and the factor 5v multiplying F3 are 



(62) 



N 
N 
N 



4sint^VK) (Q2 + Mh/)2^ 



5y = 1, + h ^ fl~ + X, 

Sy = -l^ jy + h ^ fi+ + X. 



(63) 



In principle, one can use the y dependence to determine all three of Fi,F2,Fs in a deeply 
inelastic scattering experiment. 



4.3 Space-time structure of DIS 

So far, we have used the symmetries of QCD in order to write the cross section for deeply 
inelastic scattering in terms of three structure functions, but we have not used any other 
dynamical properties of the theory. Now we turn to the question of how the scattering 
develops in space and time. 

For this purpose, we define a convenient reference frame, which is illustrated in Fig. 
Denoting components of vectors v'^ by (f f ~, v^-), we chose the frame in which 

(g+,g-,q) = i= (-g,g,0). (64) 
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Figure 19: Reference frame for the analysis of deeply inelastic scattering. 



We also demand that the transverse components of the hadron momentum be zero in our 
frame. Then 

(p^p-,p)-4^(-,^,0). (65) 

Notice that in the chosen reference frame the hadron momentum is big and the momentum 
transfer is big. 

X- X+ 



Figure 20: Interactions within a fast moving hadron. The lines represent world lines of 
quarks and gluons. The interaction points are spread out in and pushed together in x~ . 

Consider the interactions among the quarks and gluons inside a hadron, using x"*" in the 



role of "time" as in Section p.3| . For a hadron at rest, these interactions happen in a typical 
time scale Ax"*" ~ 1/m, where m ~ 300 MeV. A hadron that will participate in a deeply 
inelastic scattering event has a large momentum, p+ ~ Q, in the reference frame that we are 
using. The Lorentz transformation from the rest frame spreads out interactions by a factor 
Q/m, so that 

Ax+ ~ - X = (66) 



This is illustrated in Fig. ED 



I offer two caveats here. First, I am treating Xbj as being of order 1. To treat small Xbj 
physics, one needs to put back the factors of Xbj, and the picture changes rather dramatically. 
Second, the interactions among the quarks and gluons in a hadron at rest can take place 
on time scales Ax"*" that are much smaller than 1/m, as we discussed in Section ^ We will 
discuss this later on, but for now we start with the simplest picture. 
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Figure 21: The virtual photon meets the fast moving hadron. One of the partons is annihi- 
lated and recreated as a parton with a large minus component of momentum. This parton 
develops into a jet of particles. 

What happens when the fast moving hadron meets the virtual photon? The interaction 
with the photon carrying momentum g~ ~ Q is localized to within 

Aa;+ ~ 1/Q. (67) 

During this short time interval, the quarks and gluons in the proton are effectively free, since 
their typical interaction times are comparatively much longer. 

We thus have the following picture. At the moment x"*" of the interaction, the hadron 
effectively consists of a collection of quarks and gluons (partons) that have momenta {pf, Pi). 
We can treat the partons as being free. The are large, and it is convenient to describe 
them using momentum fractions 

ir=pt/p^, 0<e. <1. (68) 

(This is convenient because the are invariant under boosts along the z axis.) The transverse 
momenta of the partons, pj, are small compared to Q and can be neglected in the kinematics 
of the 7-parton interaction. The "on-shell" or "kinetic" minus momenta of the partons, 
p~ = p'f/{2pf), are also very small compared to Q and can be neglected in the kinematics 
of the 7-parton interaction. We can think of the partonic state as being described by a wave 
function 

^{pt,Pi;pt,P2;---), (69) 

where indices specifying spin and flavor quantum numbers have been suppressed. 

This approximate picture is represented in Feynman diagram language in Fig. |2^. The 
larger filled circle represents the hadron wave function ip. The smaller filled circle represents 
a sum of subdiagrams in which the particles have virtualities of order Q^. All of these 
interactions are effectively instantaneous on the time scale of the intra-hadron interactions 
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Figure 22: Feynman diagram for deeply inelastic scattering. 



that form the wave function. The approximate picture also leads to an intuitive formula 
that relates the observed cross section to the cross section for 7-parton scattering: 

da j^sr^ r If daai^p) 



dE' duo' 



Wf...) + (70) 



In Eq. (^), the function / is a parton distribution function: fa/h{(,, ^I'f) dC, gives proba- 
bility to find a parton with flavor a = g,u,u,d, ... in hadron h, carrying momentum fraction 
within d^ of ^ = pt /P~^- If we knew the wave functions ip, we would form / by summing over 
the number n of unobserved partons, integrating iV'nP over the momenta of the unobserved 
partons, and also integrating over the transverse momentum of the observed parton. 

The second factor in Eq. (|70D , daa/dE' du' , is the cross section for scattering the lepton 
from the parton of flavor a and momentum fraction ^. 

I have indicated a dependence on a factorization scale fip in both factors of Eq. (^). This 
dependence arises from the existence of virtual processes among the partons that take place 
on a time scale much shorter than the nominal Ax~^ ~ Q/m?. I will discuss this dependence 
in some detail shortly. 

4.4 The hard scattering cross section 

The parton distribution functions in Eq. ([70|) are derived from experiment. The hard scatter- 
ing cross sections dda^jj) / dE' doj' are calculated in perturbation theory, using diagrams like 
those shown in Fig. The diagram on the left is the lowest order diagram. The diagram 
on the right is one of several that contributes to da at order a^; in this diagram the parton 
a is a gluon. 



One can understand a lot about deeply inelastic scattering from Fig. which illustrates 
the kinematics of the lowest order diagram. Recall that in the reference frame that we are 
using, the virtual vector boson has zero transverse momentum. The incoming parton has 
momentum along the plus axis. After the scattering, the parton momentum must be on the 
cone k^j_k^ = 0, so the only possibility is that its minus momentum is non-zero and its plus 
momentum vanishes. That is 

+ g+ = 0. (71) 
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Since p+ = Q/{x^\/2) while g+ = — Q/"\/2, this imphes 

^ = x^.y (72) 

The consequence of this is that the lowest order contribution to da in Eq. ([7D| ) contains a delta 
function that sets ^ to Xbj. Thus deeply inelastic scattering at a given value of Xbj provides 
a determination of the parton distribution functions at momentum fraction equal to Xbj, 
as long as one works only to leading order. In fact, because of this close relationship, there 
is some tendency to confuse the structure functions F„(xbj, Q^) with the parton distribution 
functions fa,h{^, A^f)- I will try to keep these concepts separate: the structure functions F„ are 
something that one measures directly in deeply inelastic scattering; the parton distribution 
functions are determined rather indirectly from experiments like deeply inelastic scattering, 
using formulas that are correct only up to some finite order in a^- 



4.5 Factorization for the structure functions 

We will look at DIS in a little detail since it is so important. Our object is to derive a 
formula at lowest order in perturbation theory relating the measured structure functions for 
-\- h —>■ e~ + X via photon exchange and the parton distribution functions. 



Start with Eq. (ffOD, representing Fig. We change variables in this equation from 



{E',uj') to (xbj,!/). We relate Xbj to the momentum fraction ^ and a new variable x that is 
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just Xbj with the proton momentum replaced by the parton momentum ^p^: 



2p-q 2^p-q 



(73) 



That is, X is the parton level version of Xbj. The variable y is identical to the parton level 
version of y because p'^ appears in both the numerator and denominator: 



p ■ q C,p ■ q 
p ■ k C,p ■ k 



(74) 



Thus Eq. 



becomes 
da 



dxhj dy Jo 



/■I 1 



dxdy 



+ 0{m/Q). 



(75) 



We can calculate d&a/ {dx dy) in perturbation theory. At lowest order this is particularly 
simple, and we obtain results proportional to delta functions of Xbj/^. Using Eq. (|6^) to 
relate da/{dx]^jdy) to the structure functions Fi and F2 for 7 exchange, we obtain the simple 
lowest order results 



i^i(xbj, Q') ~ ^ E fa/h{xbi) + 0{as) + Oim/Q), 



(76) 



i^2(xbj, Q') ~ E fa/hix^i) + 0{as) + 0{m/Q). (77) 

a 

The factor 1/2 between Xbj-Fi and F2 follows from the Feynman diagrams for spin 1/2 quarks. 



4.6 fiF dependence 




Figure 25: Deeply inelastic scattering with a gluon emission. 

I have so far presented a rather simplified picture of deeply inelastic scattering in which 
the hard scattering takes place on a time scale Ax~^ ~ 1/Q, while the internal dynamics of 
the proton take place on a much longer time scale Ax~^ ~ Q/m?. What happens when one 
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actually computes Feynman diagrams and looks at what time scales contribute? Consider 
the graph shown in Fig. One finds that the transverse momenta k range from order 
m to order Q, corresponding to energy scales k~ = k^/2A;^ between k~ ~ ra^/Q and 
k~ = Q'^/Q ~ Q, or time scales Q/rn^ < Ax^ < 1/Q. 

The property of factorization for the cross section of deeply inelastic scattering, embodied 
in Eq. (|70|), is established by showing that the perturbative expansion can be rearranged so 
that the contributions from long time scales appear in the parton distribution functions, 
while the contributions from short time scales appear in the hard scattering functions. (See 
Ref. 1^ for more information.) Thus, in Fig. a gluon emission with k^ ~ is part of 



/(^), while a gluon emission with k^ ~ is part of da. 

Breaking up the cross section into factors associated with short and long time scales 
requires the introduction of a factorization scale, ^p. When calculating the diagram in 
Fig. 25, one integrates over k. Roughly speaking, one counts the contribution from k^ < 
as part of the higher order contribution to fa/h{^, f^p), convoluted with the lowest order hard 
scattering function da for deeply inelastic scattering from a quark. The contribution from 
fijp < k^ then counts as part of the higher order contribution to da convoluted with an 
uncorrected parton distribution. This is illustrated in Fig. |2^. (In real calculations, the split 
is accomplished with the aid of dimensional regularization, and is a little more subtle than 
a simple division of the integral into two parts.) 

hard scattering parton distributions 



log(l/np) log(^*) 
Figure 26: Time scales in factorization. 



A consequence of this is that both daail^p) I dE' doj' and fa/hi^yf^p) depend on ftp. Thus 
we have two scales, the factorization scale fip in ff/h{^, fJ'p) and the renormalization scale /i in 
as(/i). (When we expand da in powers of ^^(/x) then the coefficients depend on /x.) As with fi, 
the cross section does not depend on np. Thus there is an equation d{cross section) / d^p = 
that is satisfied to the accuracy of the perturbative calculation used. If you work harder and 
calculate to higher order, then the dependence on fip is less. 

Often one sets fip = fi in applied calculations. In fact, it is rather common in applications 
to deeply inelastic scattering to set /ii? = /i = Q- 



4.7 Contour graphs of scale dependence 

As an example, look at the one jet inclusive cross section in proton-antiproton collisions. 
Specifically, consider the cross section da/dEpdrj to make a collimated spray of particles, 
a jet, with transverse energy Ep and rapidity rj. (Here Ep is essentially the transverse 
momentum carried by the particles in the jet and 7] is related to the angle between the jet 
and the beam direction hj r] = ln(tan(^/2)). We will investigate this process and discuss 
the definitions in the next section. For now, all we need to know is that the theoretical 
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formula for the cross section at next-to- leading order involves the strong coupling as{fi) and 
two factors fa/h{x, ^p) representing the distribution of partons in the two incoming hadrons. 
There is a parton level hard scattering cross section that also depends on fi and fip- 




Et = 100 GeV 



Et = 500 GeV 



Figure 27: Contour plots of the one jet inclusive cross section versus the renormalization scale 
/i and the factorization scale fip- The cross section is da/dExdr] at r/ = with Et = 100 GeV 
in the first graph and Et = 500 GeV in the second. The horizontal axis in each graph 
represents Nuv = log2(2/^/-^T) and the vertical axis represents Nco = ^og2{2 ji p / Et). The 
contour lines show 5% changes in the cross section relative to the cross section at the center 
of the figures. The cm energy is ^/s = 1800 GeV. 



How does the cross section depend on fi in as{fi) and /i^ in /a//i(a;, /^f)? In Fig. I 



show contour plots of the jet cross section versus and fip at two different values of Et- The 
center of the plots corresponds to a standard choice of scales, fi = fip = Et/2. The axes are 
logarithmic, representing \og2{2fi/ Et) and \og2{2^p/ET). Thus /i and fip vary from Ep/S 
to 2Et in the plots. 

Notice that the dependence on the two scales is rather mild for the next-to-leading order 
cross section. The cross section calculated at leading order is quite sensitive to these scales, 
but most of the scale dependence found at order has been canceled by the contributions 
to the cross section. One reads from the figure that the cross section varies by roughly ±15% 
in the central region of the graphs, both for medium and large Et- Following the argument 
of Sec. this leads to a rough estimate of 15% for the theoretical error associated with 
truncating perturbation theory at next-to-leading order. 



4.8 MS definition of parton distribution functions 

The factorization property, Eq. ([70|) , of the deeply inelastic scattering cross section states that 
the cross section can be approximated as a convolution of a hard scattering cross section that 
can be calculated perturbatively and parton distribution functions fa/A{x, ^p). But what 
are the parton distribution functions? This question has some practical importance. The 
hard scattering cross section is essentially the physical cross section divided by the parton 
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distribution function, so the precise definition of the parton distribution functions leads to 
the rules for calculating the hard scattering functions. 

The definition of the parton distribution functions is to some extent a matter of conven- 
tion. The most commonly used convention is the MS definition, which arose from the theory 
of deeply inelastic scattering in the language of the "operator product expansion" [^]. Here I 
will follow the (equivalent) formulation of Ref. 0]. For a more detailed pedagogical review, 
the reader may consult Ref. [ p6[| . 



Using the MS definition, the distribution of quarks in a hadron is given as the hadron 
matrix element of certain quark field operators: 



fi/h{^,fJ'F) 



dy 

2tt 



-iip+y 



(p|^,(0,i/-,0)7+F^.(0)b). 



(78) 



Here \p) represents the state of a hadron with momentum p'^ aligned so that = 0. For 
simplicity, I take the hadron to have spin zero. The operator '?/'(0), evaluated at = 0, 
annihilates a quark in the hadron. The operator '?/'j(0, ?/~, 0) recreates the quark at x"^ = 
= and x~ = y~, where we take the appropriate Fourier transform in y~ so that the 
quark that was annihilated and recreated has momentum = The motivation for the 
definition is that this is the hadron matrix element of the appropriate number operator for 
finding a quark. 

There is one subtle point. The number operator idea corresponds to a particular gauge 
choice, A'^ = 0. If we are using any other gauge, we insert the operator 



dz-'A^{0,z~,O)ta 



(79) 



The V indicates a path ordering of the operators and color matrices along the path from 
(0, 0, 0) to (0, y^, 0). This operator is the identity operator in ^4+ = gauge and it makes 
the definition gauge invariant. 





DIS Parton distribution 

Figure 28: Deeply inelastic scattering and the parton distribution functions. 



The physics of this definition is illustrated in Fig. The first picture (from Fig. pT]) 
illustrates the amplitude for deeply inelastic scattering. The fast proton moves in the plus 
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direction. A virtual photon knocks out a quark, which emerges moving in the minus direction 
and develops into a jet of particles. The second picture illustrates the amplitude associated 
with the quark distribution function. We express F as F2F1 where 

Fi = Vexp(^-ig l^"dz-A+{0,z-,O)tay (80) 
and write the quark distribution function including a sum over intermediate states |A^): 

h/h{^,i^F) = \j% e-^^^^y-Y.{pm^.y'.^)i^F2\N){N\F,i,m\p)- (81) 



Then the amplitude depicted in the second picture in Fig. ^is {N\Fiipi{<S)\p) . The operator 
■0 annihilates a quark in the proton. The operator Fi stands in for the quark moving in the 
minus direction. The gluon field A evaluated along a lightlike line in the minus direction 
absorbs longitudinally polarized gluons from the color field of the proton, just as the real 
quark in deeply inelastic scattering can do. Thus the physics of deeply inelastic scattering is 
built into the definition of the quark distribution function, albeit in an idealized way. The 
idealization is not a problem because the hard scattering function da systematically corrects 
for the difference between real deeply inelastic scattering and the idealization. 

There is one small hitch. If you calculate any Feynman diagrams for fi/h{C,, fJ^p), you 
are likely to wind up with an ultraviolet-divergent integral. The operator product that is 
part of the definition needs renormalization. This hitch is only a small one. We simply 
agree to do all of the renormalization using the MS scheme for renormalization. It is this 
renormalization that introduces the scale into fi/hii^ l^p)- This role of np is in accord 



with Fig. roughly speaking ^p is the upper cutoff for what momenta belong with the 
parton distribution function; at the same time it is the lower cutoff for what momenta belong 
with the hard scattering function. 

What about gluons? The definition of the gluon distribution function is similar to the 
definition for quarks. We simply replace the quark field if) by suitable combinations of the 
gluon field A^, as described in Refs. [|14| and p6| . 



4.9 Evolution of the parton distributions 

Since we introduced a scale iip m. the definition of the parton distributions in order to 
define their renormalization, there is a renormalization group equation that gives the ^p 
dependence 

d d^ 

-T\ fa/h{x,fip) = ^ / -T- Pab{x/^,as{flp)) fb/h{^,fiF)- (82) 

din. jjbp ^ Jx q 

This is variously known as the evolution equation, the Altarelli-Parisi equation, and the 
DGLAP (Dokshitzer-Gribov-Lipatov-Altarelli-Parisi) equation. Note the sum over parton 
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flavor indices. The evolution of, say, an up quark 



u 



can involve a gluon {b = g) through 



the element P^g of the kernel that describes gluon splitting into uu. 



The equation is illustrated in Fig. 29 



When we change the renormalization scale fip, 
the change in the probability to find a parton with momentum fraction x and flavor a is 
proportional to the probability to find such a parton with large transverse momentum. The 
way to get this parton with large transverse momentum is for a parton carrying momen- 
tum fraction C, and much smaller transverse momentum to split into partons carrying large 
transverse momenta, including the parton that we are looking for. This splitting probability, 
integrated over the appropriate transverse momentum ranges, is the kernel Pab- 




Figure 29: The renormalization group equation for the parton distribution functions. 
The kernel P in Eq. (^) has a perturbative expansion 

P..(x/e, = pW(x/0 ^ + Pi?(x/0 f^)%.... (83) 

TT y TT J 

The first two terms are known and are typically used in numerical solutions of the equation. 
To learn more about the DGLAP equation, the reader may consult Refs. and . 

4.10 Determination and use of the parton distributions 

The MS definition giving the parton distribution in terms of operators is process independent 
- it does not refer to any particular physical process. These parton distributions then appear 
in the QCD formula for any process with one or two hadrons in the initial state. In principle, 
the parton distribution functions could be calculated by using the method of lattice QCD 
(see Ref. p^). Currently, they are determined from experiment. 



Currently the most comprehensive analyses are being done by the CTEQ ||2^ and MRS 



p7| groups. These groups perform a "global fit" to data from experiments of several different 
types. To perform such a fit one chooses a parameterization for the parton distributions at 
some standard factorization scale Hq. Certain sum rules that follow from the definition of 
the parton distribution functions are built into the parameterization. An example is the 
momentum sum rule: 

E r^ee/aM(e,/^) = l- (84) 
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Given some set of values for the parameters describing the fa/h{x,fio), one can determine 
fa/h{x,^) for all higher values of /i by using the evolution equation. Then the QCD cross 
section formulas give predictions for all of the experiments that are being used. One system- 
atically varies the parameters in fa/h{x,fJ'o) to obtain the best fit to all of the experiments. 
One source of information about these fits is the world wide web pages of Ref. [p^ . 

If the freedom available for the parton distributions is used to fit all of the world's data, 
is there any physical content to QCD? The answer is yes: there are lots of experiments, so 
this program won't work unless QCD is right. In fact, there are roughly 1400 data in the 
CTEQ fit and only about 25 parameters available to fit these data. 

5 QCD in hadron-hadron collisions 

When there is a hadron in the initial state of a scattering process, there are inevitably long 
time scales associated with the binding of the hadron, even if part of the process is a short- 
time scattering. We have seen, in the case of deeply inelastic scattering of a lepton from a 
single hadron, that the dependence on these long time scales can be factored into a parton 
distribution function. But what happens when two high energy hadrons collide? The reader 
will not be surprised to learn that we then need two parton distribution functions. 

I explore hadron-hadron collisions in this section. I begin with the definition of a con- 
venient kinematical variable, rapidity. Then I discuss, in turn, production of vector bosons 
(7*, W, and Z) and jet production. The theory for the production of heavy quarks is similar 
and I omit it. 

5.1 Kinematics: rapidity 

In describing hadron-hadron collisions, it is useful to employ a kinematic variable y that is 
called rapidity. Consider, for example, the production of a Z boson plus anything, p + p — > 
Z + X. Choose the hadron-hadron cm. frame with the z axis along the beam direction. In 
Fig. I show a drawing of the collision. The arrows represent the momenta of the two 
hadrons; in the cm. frame these momenta have equal magnitudes. We will want to describe 
the process at the parton level, a + b ^ Z + X . The two partons a and b each carry some 
share of the parent hadron's momentum, but generally these will not be equal shares. Thus 
the magnitudes of the momenta of the colliding partons will not be equal. We will have to 
boost along the z axis in order to get to the parton-parton cm. frame. For this reason, it 
is useful to use a variable that transforms simply under boosts. This is the motivation for 
using rapidity. 

Let = {q^,q~,(i) be the momentum of the Z boson. Then the rapidity of the Z is 
defined as 

y = iln(^). (85) 

The four components q^, q) of the Z boson momentum can be written in terms of four 
variables, the two components of the Z boson's transverse momentum q, its mass M, and 
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Figure 30: Collision of two hadrons containing partons producing a Z boson. The cm. 
frame of the two hadrons is normally not the cm. frame of the two partons that create the 
Z boson. 



its rapidity: 

= (ey(q2 + M2)/2, e-^^iq^ + M^)/2, q). (86) 

The utility of using rapidity as one of the variables stems from the transformation prop- 
erty of rapidity under a boost along the z axis: 

e'^q^, q^ e^'^q^ , q ^ q. (87) 

Under this transformation, 

y^y + uj. (88) 

This is as simple a transformation law as we could hope for. In fact, it is just the same as 
the transformation law for velocities in non-relativistic physics in one dimension. 




Figure 31: Definition of the polar angle 9 used in calculating the rapidity of a massless 
particle. 

Consider now the rapidity of a massless particle. Let the massless particle emerge from 
the collision with polar angle 6, as indicated in Fig. ^ A simple calculation relates the 
particle's rapidity y to 9: 

y = -\n (tan(^/2)) , (m = 0). (89) 
Another way of writing this is 

tan^ = l/sinh?/, (m = 0). (90) 
One also defines the pseudorapidity r/ of a particle, massless or not, by 

?7 = — In (tan(^/2)) or tan^ = l/sinh?]. (91) 
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The relation between rapidity and pseudorapidity is 



sinh?7 = y 1 + m^/gf. sinh?/. (92) 
Thus, if the particle isn't quite massless, 77 may still be a good approximation to y. 

5.2 7*, VF, Z production in hadron-hadron collisions 

Consider the process 

A + B ^ Z + X, (93) 

where A and B are high energy hadrons. This process and the corresponding process in 
which a W boson is produced are historically important because they are the processes by 
which the W and Z bosons were first observed ||29|| . 

Two features of this reaction are important for our discussion. First, the mass of the Z 
boson is large compared to 1 GeV, so that a process with a small time scale At ~ l/M^ 
must be involved in the production of the Z. At lowest order in the strong interactions, the 
process is q + q ^ Z. Here the quark and antiquark are constituents of the high energy 
hadrons. The second significant feature is that the Z boson does not participate in the 
strong interactions, so that our description of the observed final state can be very simple. 

In process (p3D, we allow the Z boson to have any transverse momentum q. (Typically, 
then, q will be much smaller than Mz-) Since we integrate over q and the mass of the Z 
boson is fixed, there is only one variable needed to describe the momentum of the Z boson. 
We choose to use its rapidity y, so that we are interested in the cross section da/dy. 



Pa 



Pb 




Figure 32: A Feynman diagram for Z boson production in a hadron-hadron collision. Two 
partons, carrying momentum fractions and ^b-, participate in the hard interaction. This 
particular Feynman diagram illustrates an order contribution to the hard scattering cross 
section: a gluon is emitted in the process of making the Z boson. The diagram also shows 
the decay of the Z boson into an electron and a neutrino. 

The cross section takes a factored form similar to that found for deeply inelastic scatter- 
ing. Here, however, there are two parton distribution functions: 

^ «E /' dU f dis fa/AiU /iF) h/Bi^B, /ii^) ^^^1^, (94) 

dy Th-'^A JxB dy 
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The meaning of this formula is intuitive: fa/A{^A, f^p) d^A gives the probabihty to find a par- 
ton in hadron A] fb/si^B, f^f) d^B gives the probabihty to find a parton in hadron B; daab/dy 
gives the cross section for these partons to produce the observed Z boson. The formula is 
illustrated in Fig. |3^. The hard scattering cross section can be calculated perturbatively. 



Fig. p2| illustrates one particular order as contribution to daab/dy. The integrations over 



parton momentum fractions have limits xa and x^, which are given by 



xa = 



^Jm^s, xb = e-y^Aiys. (95) 



Eq. (|9^) has corrections of order m/Mz, where m is a mass characteristic of hadronic 
systems, say 1 GeV. In addition, when daab/dy is calculated to order , then there are 
corrections of order a^~^^. 

We could equally well talk about A + B —>■ j* + X where the virtual photon decays into 
a muon pair or an electron pair that is observed and where the mass Q of the 7* is large 
compared to 1 GeV. For A + B^fi'^ + n^ + X one has the formula 

^ ^ 5 ^-^^^ ^wS^- ^''^ 

This process is historically important. Before QCD, one had partons and QED. Partons and 
QED did a good job of explaining deeply inelastic scattering. But there were other ways to 
explain deeply inelastic scattering. High mass dimuon production was investigated experi- 



mentally by Lederman et al. |30] Drell and Yan [31| proposed to explain the experimental 
results using the lowest order version of the formula above. It worked. The alternative meth- 
ods that worked for deeply inelastic scattering did not work here. This helped to establish 
the parton picture. 



5.3 Factorization is not so obvious 



The factorization formula Eq. (|9^) is supposed to hold up to m? /Q'^ corrections. This result 
is not so obvious, and in fact does not hold graph by graph. A graph for which it does not 
hold is shown in Fig. ^ Does factorization hold if one sums over graphs? The answer is 
yes, but to show this one needs to use unitarity, causahty and gauge invariance. For more 

PI- 



information, the reader is invited to consult Ref. 



5.4 Jet production 

In our study of high energy electron-positron annihilation, we discovered three things. First, 
QCD makes the qualitative prediction that particles in the final state should tend to be 
grouped in collimated sprays of hadrons called jets. The jets carry the momenta of the 
first quarks and gluons produced in the hard process. Second, certain kinds of experimental 
measurements probe the short-time physics of the hard interaction, while being insensitive 
to the long-time physics of parton splitting, soft gluon exchange, and the binding of partons 
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Figure 33: A graph for which factorization does not work. The spectator partons interact 
softly with the active partons, so that the soft part of the graph does not break up into two 
factors. 




Figure 34: Sketch of a two-jet event at a hadron colhder. The cyhnder represents the 
detector, with the beam pipe along its axis. Typical hadron-hadron collisions produce beam 
remnants, the debris from soft interactions among the partons. The particles in the beam 
remnants have small transverse momenta, as shown in the sketch. In rare events, there is 
a hard parton-parton collision, which produces jets with high transverse momenta. In the 
event shown, there are two high jets. 



into hadrons. Such measurements are called infrared safe. Third, among the infrared safe 
observables are cross sections to make jets. 

These ideas work for hadron-hadron collisions too. In such collisions, there is sometimes 



a hard parton-parton collision, which produces two or more jets, as depicted in Fig. 
Consider the cross section to make one jet plus anything else, 

A + B ^jet + X. (97) 

Let Et be the transverse energy of the jet, defined as the sum of the absolute values of the 
transverse momenta of the particles in the jet. Let y be the rapidity of the jet. Given a 
definition of exactly what it means to have a jet with transverse energy and rapidity y, 
the jet production cross section takes the familiar factored form 



One diagram that contributes to da at next-to-leading order is shown in Fig. p5| . 

What shall we choose for the definition of a jet? At a crude level, high Et jets are 
quite obvious and the precise definition hardly matters. However, if we want to make a 
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Figure 35: A Feynman diagram for jet production in hadron-hadron collisions. The leading 
order diagrams for A + B ^ jet + X occur at order a^. This particular diagram is for an 
interaction of order a^. When the emitted gluon is not soft or nearly collinear to one of 
the outgoing quarks, this diagram corresponds to a final state like that shown in the small 
sketch, with three jets emerging in addition to the beam remnants. Any of these jets can be 
the jet that is measured in the one jet inclusive cross section. 

quantitative measurement of a jet cross section to compare to next-to-leading order theory, 
then the definition does matter. There are several possibilities for a definition that is infrared 
safe. The one most used in hadron-hadron collisions is based on cones. Here I will present a 
different algorithm that is similar to the algorithms used to define jets in electron-positron 
annihilation. 



5.5 kx algorithm 

The main idea of the /ct algorithm ||32| is to modify one of the algorithms used in e^e" 



annihilation so that we use Et, rj and (j) as variables and to avoid contamination by the 
many low Et particles in the event. We choose a merging parameter R. Then we start with 
a list of "protojets" with momenta Pi, . . . as illustrated in Fig. We also start with 
an empty list of finished jets. The end result is a list of momenta pk of finished jets, ordered 
in Et- 



Figure 36: A two jet event in a proton antiproton collision. The two protojets on the lower 
left are the first to be combined. 

The algorithm can be stated very simply. See Fig. 
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1. For each pair of protojets define 



(99) 



For each protojet define 




(100) 



2. Find the smallest of all the dij and the dj. Call it d, 



3. If dmin is a dij, merge protojets i and j into a new protojet k with 



ET,k 

(f>k 




(101) 



4. If dmin is a dj, then protojet i is "not mergable." Remove it from the hst of protojets 
and add it to the hst of jets. 

5. If protojets remain, go to 1. 

Evidently, if two protojets are collincar, they will be merged right away. If one has 
vanishing momentum, it will either get merged with a protojet nearby in angle, or it will 
become a low Et jet in the final hst. Many of the jets have small Et and are really minijets, 
or just part of low Et debris. For an inclusive cross section to make n high Et jets plus 
anything else, the many low Et jets do not affect the result. For an exclusive n jet cross 
section, one would use a cutoff i?T,min- Thus in either case, low Et particles do not change 
the result. Thus the algorithm is infrared safe. 



QCD is a rich subject. The theory and the experimental evidence indicate that quarks and 
gluons interact weakly on short time and distance scales. But the net effect of these interac- 
tions extending over long time and distance scales is that the chromodynamic force is strong. 
Quarks are bound into hadrons. Outgoing partons emerge as jets of hadrons, with each jet 
composed of subjets. Thus QCD theory can be viewed as starting with simple perturbation 
theory, but it does not end there. The challenge for both theorists and experimentalists is 
to extend the range of phenomena that we can relate to the fundamental theory. 



6 Epilogue 
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